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Preface 


Under  the  impression  that  the  transmission  of  x-rays 
through  the  atmosphere  had  been  exhaustively  studied,  I  was 
suprised  to  find  that  there  is  very  little  written  in  the 
literature  which  deals  with  the  subject.  This  paper  deals 
with  the  fluence  and  transmission  of  x-rays  and  the  fluence 
or  dose  of  prompt  radiation  (neutroms  and  secondary  gamma 
rays)  and  compares  the  results  of  the  former  with  those  re¬ 
sults  predicted  by  use  of  the  Horizons  Technology,  Inc.  (HTI) 
x-ray  fluence  and  transmission  program. 

This  thesis  topic  was  one  presented  by  Dr,  Charles 
Bridgman  of  the  faculty  of  the  School  of  Engineering,  Air 
Force  Institute  of  Technology  with  the  expressed  intent  of 
validating  the  results  obtained  using  the  HTI  program  by  ob¬ 
taining  similar  results  using  mass  integral  scaling  and  the 
build-up  factor  method. 

I  wish  to  thank  Dr.  Bridgman  for  his  patience  and  heip 
in  getting  me  through  the  rough  spots  in  this  research.  I 
also  wish  to  express  my  thanks  to  my  wife,  Linda,  and  my 
children  for  the  patience  they  showed  me  even  when  I  had 
none  for  them. 


Donald  E.  Jones 
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Abstract 


This  report  is  to  validate  and  evaluate  a  Horizons 
Technology,  Inc.  (HTI)  TI-59  program  written  to  calculate 
the  free  field  x-ray  fluence  from  a  nuclear  burst.  In 
addition  to  this  validation  of  an  existing  program,  pro¬ 
grams  were  written  to  compute  mass  integral  and  prompt 
radiation  effects.  X-ray  transmission  is  calculated  using 
a  build-up  factor  method  and  is  compaied  to  results  from 
the  HTI  program.  Results  are  compared  for  black-body  temp¬ 
eratures  of  0.1,  1.0,  and  10.0  koV  and  mass  integrals  from 
10 "6  to  50,0  gm/cm2 .  The  results  compare  well  with  a  max¬ 
imum  error  of  approximately  21%.  HTI ' s  program  has  a  minor 
problem  as  the  transmission  factor  approaches  zero  so  a 
TI-59  program  is  provided  for  use  in  that  regime.  A  quick 
FORTRAN  program  is  provided  to  calculate  the  fluence  reach¬ 
ing  a  receiver.  TI-59  and  FORTRAN  programs  are  given  to 
calculate  the  mass  penetrated  and  prompt  radiation  fluence 


or  dose. 


X-RAY  FLUENCE  AND  TRANSMISSION  AND 


PROMPT  RADIATION  FLUENCE  OR  DOSE 


I.  Introduction 


Background 

The  work  reported  here  was  motivated  by  a  program 
written  by  Horizons  Technology,  Inc.  (HTI)  under  contract 
to  the  Defense  Nuclear  Agency  (contract  number  DNA  001-78* 
C-0247).  This  program,  x-ray  fluence  and  transmission 
(Ref.l),  is  one  of  a  number  of  nuclear  weapons  effects  pro¬ 
grams  developed  by  HTI.  All  of  these  programs  use  the 
Texas  Instruments  T1  -59  hand-held  programable  calculator 

(Ref. 2).  The  TI-59  stores  and  retrieves  these  programs 

t 

on/from  small  magnetic  cards. 

The  original  purpose  of  the  project  was  to  evaluate 
the  HTI  x-ray  fluence  and  transmission  program.  The  eval¬ 
uation  was  to  be  accomplished  by  comparing  the  HTI  results 
to  x-ray  transmission  factors  computed  using  a  build-up 
factor  method.  The  build-up  factors  were  calculated  for  in¬ 
finite,  homogeneous  air  using  coefficients  provided  by 
G.  Kalansky  (Ref. 3).  As  will  be  shown,  the  use  of  build-up 
factors  amounts  to  applying  the  mass  integral  scaling  approx¬ 
imation  first  described  by  Zerby  (Ref. 4)  and  recently  studied 
by  Shulstad  (Ref. 5),  The  mass  integral  scaling  approximation 
has  also  been  used  by  Murphy  (Ref. 6)  and  Eamon  (Ref. 7)  to 


calculate  free  field  neutron  and  secondary  gamma  fluences 
and/or  doses  from  a  nuclear  burst.  So,  an  extension  of  the 
project  and  the  mass  integral  scaling  technique  led  to  a 
new  program,  not  published  by  HTI ,  to  calculate  prompt  rad¬ 
iation  effects.  Since  both  the  existing  IITI  program  and  the 
build-up  x-ray  and  prompt  radiation  programs  require  the 
mass  integral  as  an  input,  a  second  new  program  was  written 
to  calculate  the  mass  integral. 

Purpose 

The  purpose  of  this  work  became  threefold.  The  first 
part  was  to  validate  the  HTI  TI-59  x-ray  flux  and  trans¬ 
mission  program  (Ref . 1 : 9- 1, 13) . 

The  second  part  of  the  purpose  was  to  write  a  new  pro¬ 
gram  to  calculate  prompt  radiation  effects.  Finally,  a  sec¬ 
ond  new  program  was  to  be  written  to  calculate  the  mass  in¬ 
tegral  . 

Method 

The  use  of  build-up  factors  to  compute  x-ray  fluence  and 
the  transmission  factor  is  developed  in  Chapter  II.  The  build¬ 
up  factor  method  (BU)  is  contrasted  to  HTI's  method  of  cal¬ 
culation,  The  programs  written  for  this  project  are  describ¬ 
ed  in  Chapter  III.  Results  of  both  methods  are  compared  in 
Chapter  IV.  Chapter  V  discusses  the  work  of  Murphy  and  Eamon, 
with  respect  to  neutrons  and  secondary  gammas,  and  the  new 
program  which  is  provided.  Chapter  VI  states  the  conclusions 


and  recommendations  of  this  work. 


Assumptions 


Two  explicit  assumptions  are  made  in  this  work.  They 


are 

1.  Kalansky's  build-up  factor  coefficients  are  applic¬ 
able  (Ref. 3). 

2.  The  concept  of  mass  integral  scaling  applies  to 
x-rays  and  prompt  radiation. 

The  assumptions  are  discussed  later  in  the  text. 
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II.  Basic  Principles 


Theory 

In  this  work,  the  HTI  program  is  evaluated  by  compar¬ 
ing  its  predictions  to  results  from  an  alternate  prediction, 
the  build-up  factor  (BU)  method.  Both  methods  use  the  same 
equation  for  fluence: 


F  « 


SfT 

'4ir'r7 


(1) 


where 

F  =  fluence  in  calories  per  square  centimeter 
S  =  source  yield  in  calories 
T  =  transmission  factor 

r  »  distance  from  source  to  receiver  in  centimeters 
f  «  x-ray  fraction  of  ahe  source  yield 
The  only  real  difference  between  the  BU  method  and  the  HTI 
model  is  in  the  calculation  of  the  transmission  factor,  T. 
For  this  reason,  a  large  portion  of  this  work  is  devoted  to 
the  calculation  of  the  transmission  ceofficient  by  each 
model , 

Transmission  Factor  for  the  BU  Method 

The  BU  method  for  calculation  of  the  transmission  fac¬ 
tor,  T  ,  :s  based  on  mass  integral  scaling  of  x-ray  trans¬ 
mission  in  a  homogeneous  atmosphere.  Mass  integral,  scaling, 
first  suggested  by  Zerby  (Ref. 4)  in  1956,  is  current !>  used 
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for  neutron  and  gamma  ray  transport  in  the  Air  Force  Wea¬ 
pons  Laboratory  computer  code  SMAUG  (Ref. 6)  and  has  re¬ 
cently  been  described  by  Shulstad  (Ref. 5).  An  explanation 
of  mass  integral  scaling  is  given  in  Appendix  A. 

This  work  made  use  of  build-up  factor  coefficients  de¬ 
termined  by  Kalansky  (Ref. 3)  to  account  for  the  increased 
number  of  x-rays  reaching  a  receiver  due  to  scattering. 
Kalansky  did  his  x-ray  transmission  calculations  using  a 
moments  method  solution  of  the  Boltzmann  transport  equation 
in  an  infinite,  homogeneous  atmosphere.  His  results  are  re¬ 
ported  in  the  form  of  BUF  as  a  function  of  mean-f ree-path 
with  x-ray  energy  as  a  parameter. 

If  the  build-up  factors  are  known  for  all  energies, 
they  can  be  used  to  compute  the  transmission  factor  for  a 
spectrum  of  x-ray  energies  by 

00 

T  =/*  P(hv)BUF(hv)exp{-[u/p]air(hv)M.I.  }d(hv)  (2) 

o 


source  to  receiver 


BUF(hv)  =  build-up  factor  for  an  x-ray  with  energy  be¬ 
tween  hv  and  hv  +  d(hv) 

Since  the  integration  of  Eq.(2)  would  be  difficult,  if  not 
impossible,  to  do  analytically,  it  is  usually  solved  by 
numerical  integration.  Eq.(2)  can  be  rewritten  in  discrete 
energy  space  as 


T 


P  BUF 
g  g 


exp[-(p/p)|ir 


M.I.] 


(3) 


where 


G  = 

P  c 
g 


(M/P) 


g 

air 


total  number  of  energy  groups 

probability  of  an  x-ray  energy  within  the  limits 
of  the  group,  g 

total  mass  attenuation  coefficient  for  air  for 
an  x-ray  energy  within  the  limits  of  the  group, 


g 

BUF  =  BUF  for  an  x-ray  of  energy  within  the  limits  of 
g 

the  group,  g 

If  the  groups  are  sufficiently  narrow,  one  does  not  have  to 
worry  about  what  value  of  m/p  and  BUF  are  appropriate  for 
each  group.  Eq.(3)  is  used  by  the  BU  method. 

Group  Probability .  The  BU  method  assumes  the  source 
to  be  a  perfect  black-body  radiator  and  that  the  source 
x-rays  can  be  represented  by  a  single  normalised  Planckian 
black-body  spectrum.  These  normalized  Planckian  functions 
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are  (Ref. 8) 


P(hv;T) 


15 

("kT)u 


(hv)3 

ehv/W. 


(4) 


where  P(hv;T)  is  the  probability  of  an  x-ray  of  energy  hv 

with  a  black-body  temperature  of  kT  keV  ,  Determination 

of  P  in  Eq.(3)  can  be  made  by  integrating  Eq.(4)  over  the 
& 

energy  range  of  each  group  (Ref, 8) 


This  integration  allows  the  computation  of  a  probability 
for  each  energy  group  being  used.  Equation  (5)  can  be  put 
in  normalized  form  by  first  defining  the  dimensionless  quan¬ 
tity,  u  ,  to  be  u=hv/kT  .  Then  Eq.(5)  can  be  rewritten  as 


(6) 


Evaluation  of  Eq.(6)  can  be  accomplished  by  integration,  read¬ 
ing  values  from  a  graph  such  as  Figure  1  (Ref. 8),  or  reading 
tables  of  the  Planck  function  and  its  integral  (Ref. 9).  In 
this  work,  the  probability  for  each  group  is  found  by  numeri¬ 
cal  integration  over  u  using  a  box  approximation  to  numer¬ 
ically  calculate  the  integral.  Note  that  due  to  the  nature 
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OCNSITY  FUNCTION,  Ff*») 


8 


CUMULATIVE  DISTRIBUTION  FUNCTION,  G(u) 


of  the  Planckian  function  and  u  ,  one  computation  of  the 
probability  for  each  group  is  all  that  is  required  since 
the  probability  of  energy  group  u^  is  the  same  regardless 
of  the  temperature  of  the  radiating  body. 

The  computation  of  the  probabilities  for  each  group 
for  the  baseline  program  is  accomplished  in  Subroutine 
PLANCK  (Appendix  C)  using  equal  Au 

Mass  Attenuation  Coefficient.  The  mass  attenuation  co¬ 
efficients  for  air  can  be  obtained  from  several  sources 
(Ref s8; 10; 11) .  The  primary  source  used  in  this  work  is 
UCRL  50174,  Compilation  of  X-ray  Cross  Sections  (Ref. 10). 

The  data  points  which  were  used  to  generate  the  polynomial 
fit  to  the  mass  attenuation  ceofficients  are  plotted  as 
circles  on  Figure  2.  The  line  in  Figure  1  was  generated 
using  the  fit  coefficients  given  in  Table  I. 

The  mass  attenuation  coefficients  were  fit  using  a 
Laurent  series  polynomial  fit  in  powers  of  1/E  ,  where 

E  is  the  x-ray  energy.  Fit  coefficients,  root-mean-square 
errors  and  plots  of  the  fits  were  generated  for  powers  of 
1/E  from  two  through  six.  The  best  fit  is  of  fifth  degree. 
The  coefficients  for  the  fifth  degree  fit  are  shown  in  Table 
I.  This  polynomial  fit  is  plotted  in  Figure  2  along  with 
the  data  points  used  to  generate  the  fit. 

The  results  of  the  mass  attenuation  coefficient  fit  are 
written  into  Subroutine  MURIIO  for  the  baseline  program  (Appen¬ 
dix  C) . 
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MASS  RTTEN  COE PP  IN  CMk*2/GM 

4  •  n*  4  * 


TABLE  X  _ _ 

Fit  Coefficients  for  Polynomial  Fit 
of  Mass  Attenuation  Coefficients 


Degree 


Coefficient 


constant 

- . 001354 

1/E 

19.7541 

-461.7632 

~/E 

1/E3 

6680.0229 

1/E4 

-3497.3643 

1/E5 

907.3575 

Mass  Integral.  The  mass  integral  is  til"  mass  ot  air 
contained  in  a  unit  area  tunnel  between  the  source  and  the 
receiver.  Mathematically,  this  is  stated  as 

r 

M.I.  =y* P(r')  dr"  (7 


where  r  is  the  slant  range  as  defined  in  Figure  3.  If  the 
air  is  assumed  to  be  exponentially  varying  in  density  accord- 
ing  to  the  local  pressure  scale  height,  then 

M-I-  =  iTS?  /  °<VeXP[-Z/HB]  l'Z'  <8) 
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r 


Figure  3.  Definition  of  Mass  Integral  Variables.  Zfi  is 
the  source  height.  ZT  is  the  receiver  height, 

0  is  the  inclination  of  the  receiver  with  re¬ 
spect  to  the  source,  and  r  is  the  slant  range, 


where 

M.I,  58  mass  integral 

6  =  as  defined  in  Figure  4 
z  =  zB  -  zT 

p(Zg)  =  density  of  air  at  the  source  height 
Hg  =  scale  height  of  atmosphere 
zB  =  source  height 
Zrp  =  receiver  height 

The  density  and  scale  height  can  be  obtained  from  I  . S .  Stan¬ 
dard  Atmosphere  1976  (Kef.  12),  Since  Zg  ,  ztJ,  ,  and  ilg 
are  all  in  units  of  kilometers  and  p  is  usually  in  units  of 
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Figure  4.  Angle  Definition. 

gm/m3 ,  a  conversion  factor  must  be  used  to  get  the  mass  in¬ 
tegral  in  units  of  gm/cm2 . 

Eq.(8)  is  directly  integrable, 

-Hn 

M. 1 .  =  ---  p(zB)[exp(-ZVHB)-l]  (9) 

Examination  of  Eq,(9)  shows  that  the  mass  integral  will  al¬ 
ways  be  positive  since  z'  is  negative  whenever  sinO  is 
negative  and  positive  when  sin6  is  positive. 

There  is  one  area  where  Eqs.(8)  and  (9)  cannot,  be  applied. 
That  is  when  the  source  and  receiver  are  co-altitude  which  im¬ 
plies  that  sinO  is  zero  and  1/sinQ  is  infinite.  For  the 
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co-altitude  case,  another  form  of  the  mass  integral  must  be 
used  (Ref. 8). 


M.I.  =  p(ZB)r 


(10) 


where  r  is  the  range  from  the  source  to  the  receiver. 

Equations  (9)  and  (10),  with  conversion  factors  to  ob¬ 
tain  the  correct  units  for  the  mass  integral,  are  written  in¬ 
to  Subroutine  MASSI  for  the  baseline  program  (Appendix  C). 

Build-up  Factors.  As  previously  stated,  the  BUF  used 
in  this  paper  are  obtained  from  Kalansky  (Ref. 3).  As 
Kalansky  points  out,  Taylor  (Ref. 13)  suggested  the  following 
equation  for  describing  build-up  factors: 


B  »  ki  exp(Ciy)  +  A2  exp(C2y) 


(11) 


where 

B  =  build-up  factor  to  account  for  the  arrival  of  scat¬ 
tered  x-rays  at  the  receiver 
y  -  number  of  mean-f ree-paths  of  source  energy 
A  2  —  1  —  A  j 

Ai,Ci,C2=  constants  determined  from  calculated  BUF 

Kalansky  determined  the  constants  A),  Cj,  and  C2  by  fitting 
his  moments  calculated  results  (Ref. 3).  The  coefficients 
Kalansky  derived  are  tabulated  in  Appendix  13  which  is  a 
table  extracted  from  Kalansky’ s  work.  When  using  this  data 
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JM ,  w  wup^ppppp . 


to  calculate  BUF,  interpolation  of  BUF,  not  Kalansky's  co¬ 
efficients,  should  be  made  when  coefficients  are  not  given 
for  the  exact  energy  of  interest. 

Final  Form.  With  the  previous  factors  in  mind,  Eq.(3) 
can  be  written  in  the  final  form  used  in  the  baseline  pro¬ 
gram 

150 

T  =  S  pe  BUFg  exp['(ll/p)Iir  M-1-]  (12) 

g=l 

Baseline  Program.  All  of  the  pieces  described  to  this 
point  are  consolidated  into  the  baseline  program  described 
and  listed  in  Appendix  C.  Additionally,  the  baseline  program 
is  set-up  to  calculate  the  4uH2  fluence.  The  program  can  be 
modified  to  calculate  the  fluence  at  a  particular  spatial 
position. 

Limitations  of  BU  Method.  There  are  two  factors  which 
must  be  considered  before  applying  the  BU  method  to  a  par¬ 
ticular  problem.  They  are: 

1.  Kalansky's  BUF  consider  mean-f ree-paths  of  up  to 

15,  His  fitted  results  did  not  converge  to  his  cal¬ 
culated  moments  results  beyond  15  mfp. 

2.  The  BUF  are  based  on  infinite  homogenous  air  and 
are  then  assumed  valid  in  exponentially  varying  air 
by  use  of  mass  integral  scaling.  This  assumption  is 
only  valid  as  long  as  the  assumption  of  mass  integral 
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scaling  is  valid.  A  literature  search  did  not  re¬ 
veal  that  any  research  has  been  done  to  determine 
the  applicability  of  mass  integral  scaling  to  the 
x-ray  transmission  problem.  Therefore,  at  this 
point,  little  can  be  said  about  the  applicability 
of  the  BU  method. 

Transmission  Factor  for  the  HTI  Model 

The  HTI  model  used  in  its  TI-59  program  is  an  empirical 

fit  to  data  extracted  from  another  source  (Ref.l),  The  HTI 

transmission  factor  equation  is 

T  -  f ( x )  B(x,6)  (13) 

and 

f(x)  =  [l  +  81.4  exp( 1 . 86x)] ~ 1  (14) 

B(x,0)  -  [A(x0-x)-l]exp[-.x(x;,-x)a]+l  (15) 

x  =  log ! o  yr  ( 16 ) 

where 

M  =  mass  integral 
9  =  black-body  temperature 

A,xo,a  =  constants  determined  from  the  empirical  lit  (Ref.l: 

9-5) 
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The  exact  origin  of  the  above  equations  is  not  clear.  The 
functions  f  and  B  are  not  defined  other  than  by  Eqs.(14) 
through  (16).  It  appears  that  f(x)  is  the  result  of  a  fit 
applied  to  a  family  of  curves  in  Ref.l  which  were  the  trans¬ 
mission  factor  excluding  build-up.  Then  B(x,0)  is  the  fit 
to  the  difference  between  f(x)  and  another  family  of  curves 
in  Ref.l  which  do  include  build-up.  Ref.l  is  not  completely 
clear  on  this  either.  If  this  surmise  is  correct,  HTI  is 
also  using  a  form  of  build-up  factor  from  an  undefined  origin 
al  source.  Thus,  we  may  be  merely  comparing  build-up  factors 
in  this  validation,  but  there  is  no  way  of  knowing  that. 

HTI's  method  allows  a  complicated  program  to  be  put  in¬ 
to  479  program  steps  for  the  TI-59,  That  was  no  small  feat 
as  can  be  seen  by  comparing  their  program  with  the  TI-59  pro¬ 
gram  written  for  this  report  and  covered  later.  However, 
there  is  a  problem  with  the  HTI  program  as  the  transmission 
coefficient  approaches  zero,  say  less  than  (0.01).  The  prob¬ 
lem  is  that  the  empirical  fit  generated  does  not  smoothly  ap¬ 
proach  zero  as  would  be  expected  but,  instead,  reaches  zero 
earlier  than  predicted  by  the  BU  method  and  then  returns  to 
significant  values  which  are  higher  than  the  values  predicted 
by  the  BU  method.  This  problem  area  will  be  discussed  furthe 
in  Chapter  IV. 
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III.  Solution  by  BU  Method 


Baseline  Program 

The  baseline  program  developed  in  this  work  uses  150 
normalized  energy  groups  based  on  .1  increments  of  u 
As  previously  pointed,  out,  the  Planckian  probability  need 
only  be  calculated  once  in  the  program  since  u  does  not 
change  with  changing  black-body  temperature.  However,  the 
average  energy  per  group  does  change.  For  the  baseline 
program,  the  average  energy  for  each  group  is  taken  as  the 
endpoint  energy  of  each  .1  increment  of  u  .  This  us¬ 
age  would  not  bo  appropriate  for  a  coarser  grouping,  but  is 
suitable  for  the  fine  grouping. 

The  baseline  program  is  written  in  the  FORTRAN  5  com¬ 
puter  programming  language  (Ref. 14).  The  algorithms  of  Ap¬ 
pendix  D  can  be  used  to  program  into  another  computer  lang¬ 
uage  if  desired. 

The  baseline  program  was  written  in  a  modular  style  to 
facilitate  changing  subroutines  (Ref, 15).  For  instance,  if 
it  is  decided  that  the  assumption  of  a  Planckian  spectrum 
for  the  source  is  inaccurate  and  another  function  better  de¬ 
scribes  what  is  happening,  then  a  subroutine  employing  the 
new  function  can  be  written  to  replace  the  subroutine,  PLi'NCK, 
currently  in  use  (Appendix  C). 
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Program  QUICK 


In  addition  to  the  150  group  baseline  program  a  pro¬ 
gram  using  ten  energy  groups  rather  than  150  groups  was  al¬ 
so  written.  This  program  was  written  to  provide  a  fast, 
easy-to-use  program  for  future  users  and  to  determine  if 
the  150  group  fine  structure  of  the  baseline  program  was 
really  necessary.  Each  group  is  constructed  to  be  of  equal 
number  density.  A  program  listing  for  program  QUICK  is  pro¬ 
vided  in  Appendix  E.  The  average  energy  for  each  group  is 
taken  as  the  energy  at  the  mid-point  probability  for  each 
group  rather  than  the  mid-point  energy  of  each  group.  In 
reality,  this  is  of  importance  only  in  the  first  and  last 
energy  groups.  In  those  cases,  using  the  mid-point  of  the 
probability  group  tends  to  weight  the  group  towards  those 
events  of  higher  probability. 

The  program  is  set  up  so  that  both  the  mass  integral 
and  the  black-body  temperature  are  read  in  as  data.  The 
program  can  be  modified  to  make  use  of  the  mass  integral 
subroutine  (MASSI)  used  in  the  baseline  program  (Appendix  C). 
Appendix  E  describes  how  the  data  is  to  be  input. 

The  results  of  program  QUICK  are  compared  to  those  of 
the  baseline  program  and  the  HTI  program  in  Chapter  IV. 

TI-59  Program 

The  TI-59  program  written  to  calculate  the  x-ray  fluence 
and  transmission  is  similar  to  the  FORTRAN  QUICK  program. 
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It  is  considerably  longer  and  more  difficult  to  use  than  the 
HTI  TI-59  program,  A  program  listing  is  given  in  Appendix  F 
as  well  as  instructions  for  its  use.  The  program  is  actually 
two  subprograms.  The  first  is  to  calculate  the  mass  integral. 
The  second  sub-program  takes  the  output  of  the  first  and 

1.  computes  the  transmission  coefficient  for  various 
groups  excluding  build-up, 

2.  then  includes  build-up  in  each  group, 

3.  then  computes  the  total  transmission  coefficient, 

4.  then  calculates  the  4irr2  fluence, 

5.  and,  finally,  calculates  the  fluence  at  a  partic¬ 
ular  spatial  distance. 

Steps  four  and  five  are  optional  depending  on  what  results 
the  user  desires. 

One  problem  with  the  sub-programs  is  their  length.  To 
execute  through  step  three  requires  826  program  steps.  With 
a  required  partitioning  of  639.39  (Ref. 2),  the  length  neces¬ 
sitates  reading  in  four  card  sides  and  executing  through 
step  two,  then  reading  in  another  card  side  and  executing 
step  three.  To  execute  the  entire  program,  excluding  the 
mass  integral  calculation,  requires  978  program  steps  and  six 
sides  of  magnetic  cards  to  be  read  in.  This  makes  the  pro¬ 
gram  somewhat  unwieldy,  but  the  results  compare  well.  See 
Chapter  IV  for  the  comparison. 


20 


IV.  Results 


Comparison  of  Results 

The  results  to  be  compared  are  those  obtained  from 

1.  HTI  TI-59  program 

2.  Baseline  FORTRAN  program 

3.  FORTRAN  program  QUICK 

4.  TI-59  QUICK  program. 

The  results  are  best  described  by  referring  to  the 
graphs  in  Figures  5,  6,  and  7.  The  plots  were  generated  us¬ 
ing  DISSPLA  (Ref. 16),  which  is  a  computer  graphics  package. 
For  each  plot,  nine  points  were  used  for  the  0.1  and  1.0 
keV  curves  ar.d  16  points  were  used  for  the  10.0  keV  curve. 
Subroutine  SPLINE  (Ref. 16)  was  used  to  smooth  the  curve. 

The  curves  all  were  generated  with  the  mass  integral  on  the 
x-axis,  the  transmission  coefficient  on  the  y-axis,  and  three 
black-body  temperatures:  0,1,  1,0,  and  10.0  keV.  Other 
black-body  temperatures  were  examined  but  not  included  to 
avoid  cluttering  the  plots.  The  other  temperatures  examined 
followed  the  pattern  shown  in  Figure  5. 

The  data  used  for  comparison  can  be  found  in  Table  II. 

Baseline  with  HTI.  Figure  5  compares  the  results  ob¬ 
tained  using  the  BU  method  with  those  obtained  from  the  HTI 
model.  As  can  be  seen,  the  results  are  comparable.  The 
largest  difference  between  the  two  is  approximately  21% 
which  occurs  at  a  mass  integral  value  of  10.0gm/cm2  and  a 
black-body  temperature  of  10.0  keV. 
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MRSS  INTEGRRI  (GM/CM**2) 


mss  INTEGRAL  (GM/CM**2) 


NASS  INTEGRAL  CGM/CN**2) 


TABLE  II 


Comparison  of  Results  for  Various  Mass  Integral 
Values  and  Black-body  Temperatures 


IE-6  a 
IE-5 
IE-4 
IE-3 
IE-2 
.1 
1. 

10. 

50 

IE-6 

IE-5 

3E-4 

IE-3 

IE-2 

.1 

1. 

10. 

50. 


Black-body 

Temperature 

(keV) 


Transmission  Factor 


IE-6 

10. 

IE-5 

10. 

IE-4 

10. 

IE-3 

10. 

IE-2 

10. 

.1 

10. 

.2 

10. 

.3 

10. 

.4 

10. 

,5 

10. 

1. 

10. 

10. 

10. 

20. 

10. 

30. 

10. 

40. 

10. 

50. 

10. 

Baseline 

HTI 

.7127 

.7824 

.4101 

.3585 

.0825 

.0483 

9E-4 

0. 

9E-10 

.0006 

0. 

.0003 

0. 

5E-5 

0. 

7E-6 

0. 

2E-7 

.9984 

.9988 

.9918 

.9926 

.9522 

.9552 

.7809 

.  7796 

.3921 

.3573 

.0633 

.0436 

IE-4 

0. 

IE-9 

.0006 

0. 

.0004 

.9998 

1. 

.9998 

1. 

.9997 

.9998 

.9988 

.9988 

.9937 

.9926 

.9750 

.9600 

.9662 

.9431 

.9610 

.9348 

.9572 

.9308 

.9541 

.9288 

.9390 

.9245 

.3346 

.4220 

.0774 

.3572 

.0188 

.0484 

.0049 

0. 

.0014 

0. 

QUICKj  TI-59 


7E-5  7E-5 


.9996 

.9956 

.9598 

.7760 

.3920 

.0606 

2E-5 


.9996 
.9956 
.9598 
.7760 
.  3920 
.0606 
2E-5 


1. 

.9998 
.9981 
.9833 
.9712 
.9621 
.  9556 
.  9505 
.9333 
.  333G 
.074  8 
.0168 
.  0038 
.0009 


.  9998 
.9981 
.9833 
.9712 
.9621 
.9556 
.9505 
.9334 
.3293 
.0727 
.0160 
.  0036 
.0008 


IE-6  =  1  x  10-6 
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QUICK  and  the  TI-59  Program  with  Baseline.  Figures  6 


and  7  compare  the  results  of  program  QUICK  FORTRAN  and  the 
TI-59  program  written  for  this  work,  respectively,  with  the 
baseline  FORTRAN  program.  The  results  are  very  close  with 
a  maximum  error  of  less  than  1%.  The  comparison  shows  that 
little  is  gained  by  using  150  energy  groups  rather  than  10 
energy  groups  as  long  as  the  groups  are  carefully  chosen. 

This  agreement  was  the  reason  QUICK  was  programmed  for  the 
TI-59  to  see  if  it  would  be  a  better  method  than  the  exist¬ 
ing  algorithm. 

Area  of  Concern 

As  mentioned  in  Chapter  II,  there  is  one  potential  prob¬ 
lem  with  the  original  HTI  program.  The  problem  occurs  as  the 
transmission  coefficient  approaches  zero  and  cannot  be  seen 
by  reference  to  Figure  5.  The  data  given  in  Table  II  shows 
the  problem,  namely  that  HTI's  transmission  factor  reaches 
zero  prematurely  and  then  jumps  back  up  to  a  value  believed 
to  be  a  little  too  high.  The  result  of  this  is  that  the  HTI 
program  can  show  the  fluence  at  a  point  to  be  zero  when  it 
actually  may  be  a  significant  value.  Also,  just  beyond  this 
mass  integral  value,  *~he  program  may  give  too  high  a  value 
for  the  fluence.  Table  III  shows  one  such  sequence  compared 
with  the  TI-59  results  of  this  work. 

Considering  a  x-ray  source  strength  of  10 12  calories 
(1  Kiloton)  and  a  black-body  temperature  of  .1  keV,  the  HTI 
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TABLE  III 


HTI  Problem  Area 
(kT  =  .1) 

Mass 

Transmission  Factor 

Integral 

HTI 

TI-59 

IE-4  a 

.0483 

.0807 

IE-3 

0. 

7E-5 

IE-2 

.0006 

0. 

.1 

.0003 

0. 

a  IE-4  =  1  x  10-4 

program  would  predict  a  4nr2  fluence  of  zero  for  a  mass  in¬ 
tegral  of  10-3gm/cm2  (receiver  roughly  0.01  kilometers  from 
the  source  at  50.0  kilometers  altitude)  while  the  BU  method 
would  predict  a  47rr2  fluence  of  6.9  x  107  calories.  Rough¬ 
ly  the  opposite  is  true  for  a  mass  integral  of  10“2gm/cm2. 

This  example  was  for  illustrative  purposes  only  since,  in 
this  case,  the  receiver  would  be  inside  the  fireball;  however, 
it  does  demonstrate  the  problem  area. 

The  problem  is  almost  certainly  the  result  of  the  empir¬ 
ical  fit  used  by  1ITI  to  fit  its  curves.  Again,  there  is  no 
problem  as  long  as  the  transmission  coefficient  is  not  close 
to  zero,  say  not  less  than  (0.01). 
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V.  Prompt  Effects 


Background 

Because  of  the  greater  number  of  interactions  between 
source  and  receiver  when  neutrons  and  secondary  gamma  rays 
are  considered,  the  method  of  build-up  factors  cannot  be  used 
for  these  radiations.  However,  the  method  of  mass  integral 
scaling  has  been  used  for  prompt  radiation  using  a  one  dimen¬ 
sional  numerical  solution  of  the  Boltzmann  transport  equation 
in  a  homogeneous  atmosphere  as  a  starting  point.  That  is, 
equation  (2)  is  replaced  by  this  solution.  The  method  used 
was  an  anisotropic  S-N  calculation  of  the  transport  of  a  sin¬ 
gle  source  neutron  with  the  energy  distribution  of  a  fission 
or  fusion  source.  The  calculation  was  done  by  Strakor  and 
Gritzner  and  first  reported  in  ORNL  4464  (Ref. 17).  The  re¬ 
sults  are  presented  as  4nr2  fluence  or  4Ttr2  dose  as  a 
function  of  mass  integral.  A  graphical  example  of  these  re¬ 
sults  is  shown  as  Figure  8.  This  figure  was  extracted  from 
Eamon's  work  (Ref. 7). 

Thus,  in  parallel  with  the  x-ray  treatment  we  have 

4irr2F  =  ST  (17) 

except  tha'c  the  4nr2  fluence  is  given  by  a  result  like  Fig¬ 
ure  8  instead  of  by  equation  (1)  through  use  of  equation  (2). 

Murphy  has  provided  a  fit  of  the  ANISN  results  of  the 
form  (Ref,6): 
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H*PI*R*»2  FLUENCE 


F(x)  =  exp(C1+C2X+C3Xt+C4X3'/2+C:;X,/2+C6X1/3) 


(18) 


where 

F(x)  =  4irr2  dose  or  fluence 
X  =  mass  integral 

=  empiric  constants  obtained  through  least 
squares  fitting  techniques 

Eamon  added  one  further  coefficient  of  the  form  C7ln  X 
inside  the  brackets  (Ref. 7)  and  provided  the  coefficients  of 
fit  to  the  neutron  and  secondary  gamma  transport  ANISN  code 
results  (Ref. 7).  The  coefficients  of  fit  for  differing 
doses  and  sources  are  in.  Table  IV. 

This  method  was  programmed  in  both  TI-59  logic  and  in 
FORTRAN  5.  The  TI-59  program  and  the  FORTRAN  5  program  ap¬ 
pear  in  Appendix  G. 

Mass  Integral  Scaling 

The  applicability  of  mass  integral  scaling  to  the  prompt 
radiation  problem  has  been  investigated  by  Shulstad  (Ref. 5). 
It  was  found  that  mass  integral  scaling  was  good  for  source 
altitudes  between  1  and  10  kilometers  (Ref. 5);  however,  at 
higher  altitudes,  the  results  obtained  could  be  as  much  as 
twice  what  is  obtained  using  Shulstad' s  two-dimensional  cal¬ 
culation.  In  a  later  study  for  the  Air  Force  Weapons  Labor¬ 
atory,  Kaman  Sciences  calculated  the  errors  generated  from 
5  to  80  kilometers  (Ref. 7). 
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TABLE  IV 


ANISf* 

1  HOMOGENEOUS 

AIR  DATA 

NEUTRONS 

DOSE 

SOURCE 

A 

B 

C 

Si  1 i con 

Thermonuclear 

- . 20795E+02  - 

. 97296E-01 

-. 17913E-04 

Tissue 

Thermonuclear 

- . 19711E+02  * 

. 98348E-01 

- . 22342E-04 

FI uence 

Thermonuclear 

-.67751E+01 

. 52690E-02 

- . 54364E-05 

D 

E 

F 

G 

.  15771E-02 

. 17924E+01  - 

. 32101E+01 

. 23746F+00 

.  18226E-02 

. 13159E+01  - 

. 151 35E+01 

- . 84022E-02 

- .  2 1468E-03 

- . 392 14E+01 

. 10875E+02 

- .  13975E  +  01 

SECONDARY  GAMMAS 

DOSE 

SOURCE 

A 

B 

C 

Si  1 i con 

Thermonuclear 

- . 25281E+02  - 

. 90163E-01 

- .  27961E-04 

Tissue 

Thermonuclear 

- . 25566E+02  - 

.  79950E-0 1 

- .  24566E-04 

FI uence 

Thermonuclear 

- . 48600E+01  - 

. 11511E+00 

-.372675-04 

D 

E 

F 

G 

.  23939E-02 

. 95659E-* 00  - 

. 11394E+01 

.  98 1 16E  +  00 

. 2 100 1 E - 02 

. 657 1 1E+00  - 

. 57599E+00 

.  9  32  7 1E  +  00 

. 31732E-02 

. 1 3350E  +  0 1  - 

. 1 30 1 1 E  +  0 1 

.95495E+00 
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TABLE  IV  Continued 


ANISN  HOMOGENEOUS  AIR  DATA 

NEUTRONS 

A  8 

C 

Silicon  Fission 

-  .  2 1780E+02  -  .  16 126E  +  00 

.  469 17F.-04 

Tissue  Fission 

-  .  18463E+02  - .  17636E  +  00 

. 53973E-04 

Fluence  Fission 

.  79627E  +  00  -.  22572E+00 

. 73701E-04 

U 

E  F 

G 

. 3886  IE-02 

. 27024E  +  0 1  -.41190E  +  01 

. 2 1249E+00 

.  44464E-02 

,  28502E  +  01  -  .  40883E+0 1 

. 17644E+00 

.  6 1 127E-02 

.  33426E  +  01  -  .  37Q18E  +  01 

. 30794E-0 1 

SECONDARY  GAMMAS 

DOSE 

A  B 

C 

Silicon  Fission 

-  .  264 16E+02  -  .  16697F  r 00 

. 56993E-04 

Tissue  Fission 

-  .26313E  +  02  -  .  16462E+00 

. 55756E-04 

Fluence  Fission 

-  .57438E  +  01  -  .  16896E+00 

. 55243E-04 

D 

E  F 

G 

. 48224E -02 

.  263GGE  +  01  -  .  35 154E  +  01 

. 109 16E  +  0 1 

.  47309E-02 

.26300E+01  - .  36007E  +  01 

.  1  1093E  +  0 1 

.47643E-02 

.  25725E  +  01  -  .  3396/E-*  01 

.  1 1089E-*  0 1 
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Results 


Results  of  the  TI-59  program  written  for  this  project 
are  compared  with  results  from  Kaman  Sciences  (Ref. 7)  in  Fig. 
9.  Plotted  on  Figure  9  is  the  4ttR2  fluence  per  source 
neutron  as  a  function  of  mass  integral.  The  circles  are  the 
points  obtained  from  Kaman  Sciences  work  (Ref. 7)  and,  in  this 
case,  can  also  be  read  from  the  TN  curve  of  Figure  8.  The 
line  is  the  fit  to  the  points  calculated  using  the  TI-59  pro¬ 
gram.  The  differences  are  negligible. 
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Figure  9.  Comparison  of  TI-59  prompt  radiation  program  results 
with  those  of  Kaman  Science  (Ref. 7). 


VI.  Conclusions  and  Recommendations 


Conclusions 

The  HTI  program  compares  well  with  the  BU  method  as 
shown  by  Figure  5  with  the  exception  of  the  problem  area 
mentioned  in  Chapter  IV. 

Additional  contributions  of  this  project  were  pro¬ 
grams  to  calculate  prompt  radiation  and  the  mass  integral. 

Such  programs  are  not  included  by  HTI  (Ref.l);  instead,  the 
mass  integral  is  left  as  an  input  variable  with  no  direction 
as  to  how  it  is  to  be  determined  and  prompt  radiation  is  not 
examined. 

Recommendat ions 

There  are  two  recommendations  to  be  made.  The  first 
applies  to  future  work.  It  is  suggested  that  a  study  be 
undertaken  to  determine  the  validity  of  using  mass  integral 
scaling  in  the  x-ray  problem. 

The  second  recommendation  applies  to  the  use  of  the  HTI 
TI-59  program  versus  the  use  of  the  program  written  for  this 
work.  Due  to  the  HTI  program's  ease  of  use,  combined  with  re¬ 
sults  which  are  comparable,  it  is  recommended  that  the  HTI 
program  be  used  to  calculate  the  transmission  coefl Lcient  and 
fluence  as  long  as  the  transmission  coefficient  is  not  less 
than  approximately  (0.01).  If  the  transmission  coefficient  is 
less  than  0.01,  then  it  is  felt  that  the  TI-59  program  written 
for  this  work  should  be  used  to  avoid  the  curve  fitting  prob¬ 
lem. 
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Appendix  A:  Mass  Integral  Scaling 


The  mass  scaling  law  can  be  stated  in  the  following 
manner  (Ref. 5):  In  an  infinite  homogeneous  medium  with  an 
isotropic  point  source,  the  4iTr2  fluence  is  a  function 
only  of  the  density  and  range  from  source  to  receiver: 

MI  =  pR  .  An  assumption  must  be  made  to  allow  the  appli¬ 
cation  of  the  mass  scaling  law  to  the  actual  case.  The  as¬ 
sumption  made  is  that  if  the  mass  range  between  a  source  and 
a  receiver  in  infinite  homogeneous  air  is  equal  to  the  mass 
range  in  variable  density  air,  then  the  4irr2  fluences  will 
be  equal.  That  is,  points  from  the  infinite  homogeneous  air 
case  will  map  directly  to  points  in  the  variable  density  air 
case  at  identical  mass  integral  values.  For  a  derivation  of 
the  mass  scaling  law,  see  Shulstad  (Ref .5:86-89). 

While  Shulstad  (Ref. 5)  has  examined  mass  integral  scal¬ 
ing  as  it  applies  to  neutrons  and  gamma  rays,  no  references 
were  found  in  the  literature  to  the  same  type  work  being  ac¬ 
complished  for  x-rays. 
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Appendix  B:  Kalansky's  Coefficients  (Ref .3) 


I 


I 


Energy  in  I;eV 

A1 

12 

-0.227 

14 

-0.370 

16 

-0.323 

18 

-0.634 

20  . 

-1.072 

22 

-1.048 

'  24 

-1.740 

26 

-2.673 

28 

-2*664 

>0 

-6.038 

32 

-8.805 

34 

-0.504 

36 

-75.83 

38 

-20.03 

40 

-16.94 

45 

.  ‘  14.59 

50 

11.31 

55 

,  109.2 

60 

-11.05 

70 

-114.1 

80 

-113.1 

100 

-10,93 

120 

-8.153 

150 

13*14 

200 

-88.92 

250 

-6.308 

300 

19.89 

350 

72.37 

.  400. 

-6.063 

.  ■  500 

16.29 

600 

-57.53 

750 

-17.20 

A2 

Ci 

C2 

1.227 

-0.400 

0.000 

1070 

-0.400 

0.000 

1023 

-0.680 

0.020 

1*634 

-0.460 

0.020 

2.072 

-O.36O 

0.020 

2.048 

-0.480 

0.040 

2.748 

-0.340 

0.040 

3*673 

-0.260 

0.040 

3.664 

-0.300 

0.060 

7.038 

-0.140 

C.C40 

9,805 

-0.100 

0.040 

9.504 

-0.100 

0.060 

76.83 

0.000 

0.020 

21.03 

-0.02c 

0.C60 

17.94 

-0.020 

0.060 

-13.59 

0.120 

-0.020 

-10. 31 

0.160 

-0.044 

-108.2 

0.120 

0.100 

12.05 

0.000 

0.200 

115.1 

0.140 

0.160 

114.1 

0.160 

0.180 

11.93 

0.060 

0.260 

9.153 

0.020 

0.280 

-12.14 

0.260 

0.100 

89.92 

0.180 

0.200 

7008 

0.000 

0.260 

-18.89 

C.2C0 

0.120 

-7107 

0.160 

0.140 

7.063 

0.  coo 

0.220 

-15*29 

0.160 

0.080 

57.58 

0.100 

0.120 

lo.20 

O.COG 

0.120 
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Appendix  C:  Baseline  Program  Listing  aid  Use 


The  baseline  program,  with  all  its  subroutines,  is  list¬ 
ed  on  the  next  few  pages.  The  program  was  written  with  future 
users  in  mind;  that  is,  comments  are  provided  as  documentation 
to  make  the  program  easier  to  use.  Additionally,  Table  V  pro¬ 
vides  a  list  of  the  input  variables  and  their  units. 


TABLE  V 

Variable 

Representing 

Units 

K,  AX|  A2,Cj,C2 

Read  in  of  Kalansky's  data 
(Appendix  B) 

K  in  keV 

N,  NN,  L 

Dummy  variables  (see  comments 
in  program) 

— 

YLD 

Yield  of  source 

Kilotons 

XF 

X-ray  fraction  of  source 

— 

XKT 

Black-body  Temperature 

keV 

RHO 

Density  of  air  at  source 
altitude 

gm/m 

R 

Distance  from  source  to 
receiver 

kilometers 

ZB 

Source  height 

kilometers 

ZT 

Receiver  height 

kilometers 

1IB 

Scale  height  of  atmosphere 

kilometers 

ANGLE 

As  defined  in.  Figure  4 

degrees 

CM 

Slant  range  from  source  to 
receiver 

cent imetors 

- 
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PROGRAM  TRAMS 

D I MENS  I CN  F G (Oil  SO } ,  PC ( 0 1  150 ) 

DIMENSION  HNU( 150) , 2(01 150) .Y< 150) »X(150) ,XMU< 150) 

D I  MENS ION  TO ( 0 : 1 50 > , B23 ( 1 50  5 

DIMENS  ION  3 ( 1 50 ) » TY ( 0 1  1 50  ) , YY  <  0 1 1 50  > , YYY ( 150 ) » SB ( 0 1  1 50 ) 
DIMENSION  V(0: 150 . A(0: 150) . AAA ( 150) 

DIMENSION  K ( 32 ) » A1 (32 ) , A2 ( 32 ) » Cl ( 32 ) » C2 (32) 

READ  IN  BUILD-UP  FACTORS  FOR  LATER  USE. 

SOURCE  USED  IN  THIS  WORK  I  UNPUBLISHED 
MASTERS  THESIS,  'X-RAY  BUILD-UP  FACTORS', 

KALANSKY.  G.M.  AFIT,  SCHOOL  OF  ENG, 

DEC  1S7S. 

DO  2  1=1,32 

READ  (■«•»#,  END  =  S9S  )  K  (  I )  ,  A1  ( I )  ,  A2  <  I  )  ,  C 1  ( I )  ,  C2  ( I ) 
CONTINUE 

CA;_L  PLANCK  (PG,  PC) 

KK  *  0 
CONTINUE 
KK  =  KK  +  1 

N  IS  USED  TO  TELL  THE  PROGRAM  HOW  THE 
MASS  INTEGRAL  IS  TC  BE  HANDLED.  N*1 
IMPLIES  THAT  THE  TARGET  AND  BURET  ARE 
CO-ALTITUDE.  iM  =  2  IMPLIES  T.-H  MASS  INT¬ 
EGRAL  VALLES  ARE  TO  BE  READ  IN'.  ANY 
OTHER  VALLE  FOR  N  IMPLIES  THE  M.I.  IS  TC 
BE  COMPUTED  AND  THE  TARGET  AND  BURST  ARE 
NOT  CO-ALTITUDE. 

NN  IS  USED  TO  TELL  THE  PROGRAM  WHE"HER 
OR  NOT  TO  CALCULATE  THE  4  PI  R**2  FLUENCE 
(NN= 1  IMPLIES  CALCULATE). 

IF  1=0,  THEN  THE  PROGRAM  WILL  NOT  IN¬ 
CLUDE  BUILD-UP  IN  THE  TRANSMISSION  CO¬ 
EFFICIENT  CALCULATION. 

READ  (*»■», END  =  S9S)N 
READ < * •  *» END = 3S2 )  NN 
READ (*,*»END=E£S)  L 
CALL  MASS  I  (,N  ,  XM I ) 

IF(NN.EO.l)  THEN 

READ  v  *  »  * , £nD  =  SSG )  YLD,XF 
END  IF 


THIS  PORTION  COMPUTES  A  MASS  AT“ENUAT I CN  CGI 
EVERY  ENERGY  CROUP.  U=HNU/KT. 

THIS  PROGRAM  IS  SET  UP  7C  USE 
SEVERAL  VAi_l.ES  FOR  ELACR  EC  DR*  TZ'IP. 

THE  PROGRAM  CAN  EE  MODIFIED  TO  READ 
IN  BLACK  BODY  TEMPS  OF  INTEREST. 
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do  "co  L_=:,iz 

a  r  t  w w  i  i  *  i  •'  r*« 

XKT  =  .1 

ELSE  I F  <  LL . GT . 1 . AND . LL . LE . 1 0 )  THEN 
XKT  =  XKT  f  . l 
ELSE  IF(LL.GT.IO)  THEN 
XKT  *  LL  -  9. 

END  IF 

.PRINT*,  'BLACK  BODY  TEMP*  ',XKT.. 

PRINT*,  '  ' 

U=.l 

DO  3  3  *1,150 
HNU ( I >  =U*XKT 
U=U+. 1 
3  CONTINUE 

CALL  MURHO < HNU , XMU ) 

CALL  MFP  <  XMI , XMU , X ) 

IF(L.NE.O)  THEN 

CALL  BUr  ( HNU ,  3  >  IE0F,X,K,AI  » A2 »  Cl »  C2 , KK  »  BB  »  BBB, PG ) 

IFdEOF.EG.O)  GO  TO  SS9 

ENDIF 

CALL  TRNS ( TY , YY . TG , Z , V , X , Y , PG , 3 , YYY , L ) 
lF(NN.cG.l)  THEN 

CALL  FLUENCE i TO ( 1 50 ) , VLB, XF > 

ENDIF 

PRINT*,  '  ' 

PRINT*,'  ' 

300  CONTINUE 
GO  TO  10 

338  PRINT*, 'ERROR!  END  OF  DATA  AT  IMPROPER  TIME' 

SS9  STOP 
END 


SUBROUTINE  PLANCK < PG , PC ) 

DIMENSION  PGtOIlSO),  PC<0:i50> 

rc(0)=o. 

A=15./( (ACOSC-l . ) )**4> 

PC(C)=0. 

U=.  1 

DO  1  1=1,150 

PG< I >= A* ( U**3 )/( EXP ( U ) - 1  ) 

PC ( I ) =  < <PG< I )+PG( 1-1 ) >/2. >*. 1 
U=U+.l 

1  CONTINUE 
RETURN 
END 
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IF  N  =  1 ,  READ  IN  THE  D l N 3 1 T V  OF  AIR  AT 
BURST  HEIGHT  IN  GM/M**3  AND  THE  DISTANCE 
FRON  BURST  TO  TARGET  IN  KILOMETERS. 

IF  N  =  2,  READ  IN  THE  MASS  INTEGRAL  VALUES 
IN  GM/CM**2. 

IF  N  NOT  EQUAL  TC  1  OR  2,  THEN  VALUES 
MUST  BE  READ  IN  (IN  ORDER)  FOR: 


ZB*HGT  OF  BURST  (KM) 

ZT=HGT  OF  TARGET  (KM) 

RHO=DENSITY  AT  Z3  (GM/M**3) 

HB=SCALE  HGT  OF  ATMOSPHERE  (KM) 

( SOURCE lU.S.  STANDARD  ATMOSPHERE) 

ANGLE* ANGLE  BETWEEN  HORIZONTAL  AND 

STRAIGHT  LINE  FROM  BURST  TO  TGT 
(DEGREES).  ANGLE  IS  MINUS  IF 
THE  TGT  IS  LCaER  IN  ALTITUDE 
THAN  THE  BURST. 

IF(N.EQ.Z)  THEN 

READ  ( *  ,  #  , EMD  =  2S8 )  XMII 
XMI  =  XMII 

PRINT# , 'THE  MASS  INTEGRAL  WAS  AN  INPUT  VALUE' 
GO  TO  50 
ENDIF 

IF(N.EG. 1 )  GC- TO  40 

READ ( *  ,  *  , END  =  99S )  ZB  fZT  ,RHO  ,HB  f ANGLE 
PRINT*, 'BURST  ALT  =  ',ZE,'  TGT  ALT  *  ' , ZT , 

*'  DENSITY*  ',RHO,'  SCALE  HGT  AT  BURST  =  ' , HB 
PRINT*, 'ANGLE  OF  TGT  FROM  BURST  =  ', ANGLE 
ANGLE  =  ANGLE  *  .01745 
GO  TO  41 

40  READ ( * , * , END=93B )  RHO.R 
XMI *RKO*R/ 1 0 . 

GO  TO  50 

41  CONTINUE 
Z1=0. 

Z2=ZT-ZB 

XMI * ( -HB ) * ( RHO ) *  < 1 000 , ) * ( 1 . E-4 ) * ( EXP ( -Z2/H5 5  - 1 . )/ 
*SIN( ANGLE ) 

50  PRINT#, 'THE  MASS  INTEGRAL  IS  ',XMI,'  GM/CM*^C' 

GO  TO  51 

9S3  STOP  'END  OF  DATA  IN  MASSI  AT  WRONG  TIME' 

51  RETURN 
END 
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w  L.  "•  \  G  w  I  i  N  C  2  >_•  r  (t'.'toiOi  i.  C.  U  •*  r  /  \  f  A  >  1-1  -  »  n  _  »  U  *  ;  I - »  A  > 

*  3  3  •  £•  L  8  ,  r  G  ) 

DIMENSION  3(150)  ,HNU(  150>  ,  33(01  150)  ,E33(  150) 
DIMENSION  K 02) ,Ai (32) . A2<22) »C1 (32) ,02(02) 
DIMENSION  }((  1?0)  »PG(0:  lfO) 

IF ( HNU (150) .LT.12)  THEN 
DO  31  1=1,150 
31  B ( I )  =  1. 

GO  TO  35 
ENDIF 
N  =  12 

DO  1  1=1,150 

IFtHNUC)  .LT.  12)  GO  TO  12 
IF(HNU(I) .GT.750)  GO  TO  48 
DO  5  J= 1 ,32 

IF ( K ( J ) . GE . HNU ( I > )  GO  TO  3 
IF ( K ( J ) , LT.N)  GO  TO  5 
All  ■  A 1 ( J  > 

A22  =  A2 (J) 

C 11  ■  Cl ( J ) 

C22  =  C2 ( J ) 

N  ■  K(J) 

5  CONTINUE 

3  A  =  K  (  J  )  -  N 

DD  =  HNU ( I )  -  N 
D  *  DD  /  A 

BH  «  (A1(J)*EXP(X(I)*C1(J)))+  ( A2 ( J  >  * 

*  EXP ( X( I >*C2 ( J )  )  ) 

BL  *  (All  *EXP  ( X ( I ) *C1 1 )  >  +  (A22*EXPO<(  I  )*C22)  > 
B( I )  *  iBh  -  3L)*D  +  BL 
GO  TO  10 
12  B(I)  «  1. 

10  CONTINUE 

1  CONTINUE 
35  CONTINUE 
BUFF  =  0. 

—  BB(O)  =  -0. 

DO  600  1=1,150 

BB ( I )  ■  B ( I )  *  P  G (  I ) 

BBB(I)  =  ((33(1 )+BB( I-l))/2.)*.l 
BUFF  =  BBB(I)  *  BUFF 
BOO  CONTINUE 

PRINT*, 'TOTAL  BUF  =  ' , BUFF 

GO  TO  50 
49  IEOF  =  0 

PRINT*,  'ERROR!  END  OF  DATA  AT  IMPROPER 
♦TIME  IN  BUF  ' 

GO  TO  50 

48  PRINT*, 'HNU  IS  TOO  LARGE  F JR  KALANSKY3 

*  BUILD-UP  >  750' 

50  RETURN 

END 
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*?G»3»YYY.L) 

D I  MENS I ON  TV ( 0 : 1 50 ) , YY ( 0 : 1 00 > , TG ( u : 1 50  > , 

*2  <  o : i so ) . v  <  o : : oo  > , x ( i so  > ,  y  d so > 

*  f  PG ( 0 : 1 50 ) »  3 ( 150 ) , YYY ( i 00 ! 

TY(O)  =  0. 

YY ( 0 )  =  0. 

TGCO) -0 . 

Z(0)=0. 

V(0)  =  0. 

DO  5  1=1,150 

Y<I)*EXP(-X< I) ) 

IF ( Y( I ) . LT . 1 . E~20 )  Y ( I ) =0 . 

YY ( I )  =  PG  (I  )  *  Y  < I ) 

YYY ( I )  =  ( ( YY ( I )  +  YY(I-l))/2.)*.l 
TY ( I )  =  YYY ( I )  +  TY(I-l) 

IF<L . EG . 0 )  THEN 

Z(I)  =  Y(I)  *  PG < I ) 

ELSE 

Z  < I )  =  Y(X>»3(X)*PG(X) 

END  IF 

v<X)  =  (C(i)+:(M))/:.)*.i 

TG(I)  =  0(1)  +  TG(I-l) 

5  CONTINUE 

IF<L.EC5.0)  THEN 

PRINT*?  'DIRECT  TRANSMISSION  COEFFICIENTS 

*  TG( 150) 

ELSE 

PRINT*, 'TRANSMISSION  COEFFICIENT  =  ',TG(150> 
ENDIF 

PRINT*, 'INTEGRATED  EXP  ATTEN=  ',TY(150) 

RETURN 

END 


SUBROUTINE  MURHO ( HNU , XMU ) 

DIMENSION  HNU( 150) ,XMU( 150) 

DO  1  1=1,150 
E  -  1.  /  HNU ( I ) 

XMU(I)  =  - . 00 14+ ( iS.7541*E)-(4GI .7332+E**2) * 

*  ( G E 8 0 .  0222*E**3  )  -  (  34  S7  .  j343^E**4  ) 

*  + ( 907. 3575 *E**5 ) 
l  CONTINUE 

RETURN 

END 


45 


nooonnon 


CGMPUTES  KEAN  FREE  PATH  FOR  VARYING 
X-RAY  ENERGIES . 

IF  X  IS  >  100,  THE  NEGATIVE  EXPONENTIAL 
OF  X  RESULTS  IN  AN  UNDERFLOW.  THEREFORE, 
X  IS  JUST  SET  TO  100  IF  IT  EXCEEDS  100. 

DO  B  1=1,150 

X(I)=XMU(I)*XMI 
IF(X< I ) .GT. 100)  X(I)=100. 

6  CONTINUE 
RETURN 
END 


SUBROUT  I NE  FLUENC E ( TRANS , YLD , XF ) 

XYLD  =  YLD  *  XF 

PRINT*,  'X-RAY  YIELD  IN  KILOTONS=  ' , XYLD 

YLDD  =  XYLD  *  1.E12 

FL  *  YLDD  *  TRANS 

PRINT*, '4  PI  R**2  FLUENC£=  ' , r L 

RETURN 

END 
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Appendix  D:  Algorithm 


The  intent  of  this  appendix  is  to  provide  the  user  with 
an  algorithm  which  can  be  programmed  into  any  suitable  com¬ 
puter  language.  The  algorithm  is  based  on  the  numerical  in¬ 
tegration  of  Eq.(12)  with  the  involved  factors  properly  cal¬ 
culated.  Proceed  as  follows: 

1.  Decide  whether  mass  integral  values  are  to  be  input 
or  calculated.  If  input,  proceed  to  step  two.  If 
calculated,  use  either  Eq.(9)  or  (10)  as  the  case 
dictates.  Density  and  scale  height  can  be  found  in 
Ref. 12  and  source  height,  receiver  height,  and  theta 
or  range  should  be  provided  by  the  problem  definition. 

2.  Determine  the  probability  of  each  group  using 

Pg  -  (IS/oj^l]  (19) 

Store  the  results.  For  the  baseline  program,  u  was 
between  0  and  15,  inclusive,  and  was  incremented  by 
0.1  for  each  group. 

3.  Input  the  Planckian  black-body  temperature.  Deter¬ 
mine  the  energy  of  each  group  by 

(hv)  =  u  *  kT  (20) 

o  o 

Store  for  future  use. 

4.  Calculate  the  mass  attenuation  coei'ficii  iu  !or  each 


47 


group  using 


E  =  1/ (hv )  (21) 

b  b 

(y/p)  =  -.00 14 +(19. 754 1*E  )- (461 . 7632*E  2) 

6  6  6 

+  (6680. 0229*E  3 )- ( 3497 . 3643*E  4 ) 

6  6 

+(907.3575*Egs)  (22) 

Store  each  as  it  is  calculated. 

5.  Calculate  the  mean-f ree-path  for  each  group. 

(MFP)g  «-  (y/p)  *M.I.  (23) 

Store  the  results. 

6.  Input  Kalansky’s  build-up  coefficients.  Calculate 

the  build-up  factor  for  each  group. 

a.  If  (hv)  is  less  than  12  keV,  then  set 

S 

(BUF)  equal  to  1.0 
6 

b.  If  (hv)  is  greater  than  750  keV,  stop. 
Kalansky's  coefficients  are  not  valid  above 
750  keV. 

c.  If  (hv)  is  greater  than  12  keV  and  less 

S 

than  750  keV,  interpolate  for  (BUF)  be¬ 
ta 

tween  given  energies  for  which  Kalansky  pro¬ 
vides  coefficients.  Note;;  Interpolate  be¬ 
tween  BUF  not  coefficients.  As  used  below, 
h  =  high  and  1  =  low 
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(BUF)h  : 
and 

(BUF),  « 

(hv)T  « 

(hv)p  - 

(BUF)g  = 

7. 


(Ai)h  exp 


'(C,)h(MFP) 


g 


+(A2)h  exp 


1( 


C2)h(»'FP)g 


(24) 


(Ai)i  exp£(C,),(MFP)gJ+(A?)i  exp  ^(C2  )  ,  (MFP)  1 


(25) 


s  <*»h  -  <hv>i 


(26) 


=  (hv)  -  (hv)4 


(27) 


[(BUF)h-(BUF)J*^hv)p/(hv)J+(BUF)Jl  (28) 


Store  the  results. 

Calculate  the  transmission  factor,  first  for  each 
group,  then  sum.  This  step  makes  use  of  the  pre¬ 
vious  six  steps. 


(29) 


(Trans)  =  P  *  Y  *  (BUF) 
g  B  g  g 


(30) 


Finally,  integrate  and  sum.  The  baseline  program 
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uses  a  box  approximation  to  evaluate  the  integral. 


Trans 


G 


Trans ) 


g+(Trans) 


(31) 


8.  'The  4ur2  fluence  is  found  by  multiplying  Eq.(31) 

by  the  x-ray  fource  strength  in  calories,  XS. 

4ttt2  fluence  =  XS*Trans  (32) 

9.  Finally,  the  fluence  is  found  by 


Fluence  = 


XS*Trans 

4irr2 


(33) 


where  r  is  the  distance  from  the  source  to  the 


receiver  of  interest  in  centimeters. 


Appendix  E:  Program  QUICK  Listing 


As  previously  explained,  program  QUICK  is  a  FORTRAN 
program  which  uses  10  equi-probability  groups  rather  than 
150  equal  energy  spacing  groups.  Program  QUICK  is  designed 
to  be  quick  running  and  simple.  A  program  listing  is  pro¬ 
vided  on  the  next  two  pages.  For  convenience,  Kalansky's  co 
efficients  are  read  in  as  the  first  piece  of  data  in  the  pro 
gram  since  an  interactive  computer  system  was  used  to  facili 
tate  program  modification.  The  coefficients  could  be  put  in 
DATA  (Ref. 14)  statements  if  so  desired. 

All  factors  should  be  easily  identifiable  with  the  pos¬ 
sible  exceptions  of  A,  XMI ,  XKT,  YLD,  XF,  and  CM.  They  are 
A  -  the  mid-probability  value  of  u  for  each  group 
XMI  -  the  mass  integral 

XKT  -  the  Planckian  black-body  temperature 
YLD  -  yield  in  kilotons 
XF  -  x-ray  fraction  of  YLD 

CM  -  range  from  source  to  receiver  in  centimeters 
Program  QUICK  is  written  with  the  mass  integral  as  an 
input;  however,  it  could  be  modified  to  compute  the  mass;  in¬ 
tegral  using  subrout ine  MASSI  (Appendix  C). 
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M^**4»A4iiAJ>**OOUUUUUUUUNMNHNNNNMM^i-mHmmmh>mm 

OU(S>iOHn>UN^O(OOB>i(MQ»UN^OtOCO>ltnUI«UrO^O(90>l(nCI^UN^Ot30lNl<nU^UN 


1 


2 


3 


5 

6 


11 

9 
4 

10 


DIMENSION;  A< 10) ,K<32) ,A1 (32) , A2<32) ,Ci <32) 
DIMENSION  C2 ( 32 )  ,  HNU  (  10 )  ,XMU( 10)  ,X(  10)  ,3(10) 
D I  MENS I  ON  Y ( 0 1 1 0 ) ,70(01 10) 

DO  1  1=1,32 


READ  (•*»*,  cND  =  9S9  )K  ( I )  ,  A1  ( I )  ,  A2  ( I )  ,  Cl  ( I )  ,  02  ( I ) 
CONTINUE 

DATA  A/1. 1 ,1 .8,2.35,2.8,3.25,3.75,4.3,5. ,5. 2,7.7/ 


READ  <  * ,  #  ,  £ND  =  9SS )  XMI ,  XKT  ,  YLD  ,  )(F  ,  CM 
DO  2  1*1,10 

HNU(I)  =  A < I )  *  XKT 
E  =  1.  /  HNU(I) 

XMU(I)  =  -.0014  +  (19.7541*E>  -  <461.7632* 

*  E*E )  +  <GB80.0222'f-E->*3)  -  (3497. 3G43 

*  *E**4 )  +  (907.3575*E^*5) 

X(I)  =  XMU(I)  *  XMI 

IF(X(I)  .GT.100)  >(<I)  =  100. 

Y<  I )  =  EXP(-X(D) 

CONTINUE 

IF <HNU( 10) .LT. 12)  THEN 
DO  3  1=1,10 
B ( I )  «  1. 

GO  TO  10 
ENDIF 
N  =  12 
DO  4  1=1,10 

IF ( HNU ( I ) . LT . 1 2 )  GO  TO  11 
IF < HNU < I ) . GT . 750 )  GO  TO  99 
DO  5  J= 1 , 32 

IF < K  < J) .GE.HNUt I ) >  GO  TO  G 
IF(K< J) .LT.N)  GO  TO  5 
All  =  Al ( J ) 

A22  =  h2 ( J ) 

Cll  =  C1(J) 

C22  =  C2 ( J ) 

N  =  K  <  J  > 

CONTINUE 
.  AA  =  K ( J )  -  N 
DD  =  HNU ( 1 )  ~  N 
D  =  DD  /  AA 

BH  =  (A1<J)*EXP(X(!)*C1(J)) )  +  ( A2 ( J ) 

*.  *EXP(X(I)*C2(  J)  )  ) 

BL  =  <A11*EXP(X(I )*C11 > )  +  (A22«EXP(X( I ) 

*  *G22 ) ) 

B ( I )  =  (BH  -  BLMD  +  BL 
GO  TO  9 
B  ( I )  =  1. 

CONTINUE 

CONTINUE 

CONTINUE 
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51  Y(0)  =  0. 

52  TG ( 0 )  =  0. 

53  DO  7  1=1,10 

54  IF(Y(1) .LT.l.E-20)  Y < I )  =  0. 

55  Yd)  »  Yd)  *  .1  *  B<I> 

56  TG< I )  =  Y(i)  +  TG(I-l) 

57  7  CONTINUE 


58  XYLD  =  YLD  *  XF 

58  YYLD  =  XYLD  *  1.E12 

60  FL  =  YYLD  *  TG(10)  / (4. *AC0S<-1 . ) *< CM*CM) > 

61  PRINT*, 'FOR  A  BLACK  BODY  TEMPERATURE  OF  ',XKT,'  KEV, 

62  *  A  MASS  PENETRATED  OF  ',XMI,'  GM/CM**2,  A  WEAPON 

63  *  YIELD  OF  ' , YLD ,  '  KILOTGNS,  A  X-RAY  FRACTION  OF  ' 

64  .  #,XF»'  AND  A  DISTANCE  FROM  BURST  TG  TARGET  GF  ', 

65  *CM»  *  CM,  THE  TRANSMISSION  COEFFICIENT  IS  #,T0<10> 

66  AND  THE  FLUENCE  IS  '  , r L 

67  GO  TO  1 

69  S3  PRINT*, 'KNU  IS  TOO  LARGE  FOR  KALANSKYS  BUF-  >  750' 

69  993  STOP 

70  END 
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Appendix  F:  TI-59  Programs 


Comments 

The  program  listings  provided  on  the  next  nine  pages  are 
those  written  during  this  research  to  calculate  the  trans¬ 
mission  factor,  47rr2  fluence,  and  fluence  at  a  particular 
point  using  the  TI-59.  For  a  recommendation  on  when  to  use 
this  group  of  programs,  see  Chapter  VI. 

The  program  listed  on  page  55  is  to  calculate  the  mass 
integral.  Such  a  program  is  not  provided  in  HTI's  work.  The 
program  on  pages  56  through  60  is  to  calculate  the  transmission 
factor  for  each  group.  The  program  is  not  complete  by  itself. 
It  needs  at  least  a  portion  of  the  program  (the  first  186 
steps)  on  pages  61  through  63  to  determine  the  transmission 
factor  for  the  situation  being  investigated.  The  other  por¬ 
tions  of  the  program  listed  on  pages  61  through  63  are  to  com¬ 
pute  the  4irr2  fluence  and  the  fluence  at  a  particular  point. 

Because  the  use  of  the  TI-59  programs  provided  in  this 
work  is  somewhat  cumbersome,  instructions  for  their  use  are 
provided. 

Instructions  for  Use 

Note:  Partitioning  is  639.39  for  all  TI-59  programs  writ¬ 
ten  for  this  work  (Ref,2). 

1.  If  mass  integrals  are  already  known,  skip  to  step 
two.  If  not,  read  in  the  one  card  side 


Mass  Integral  Program 
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Transmission  Factor: 
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for  the  mass  integral  program  (page  55).  Then 

a.  If  the  source  and  receiver  are  co-altitude : 
Press  A 

Enter  density  in  gm/m3  Press  R/S 
Enter  range  in  km  Press  R/S 

*  Read  mass  integral 

b.  If  the  source  and  receiver  are  not  co-alti¬ 
tude: 

Press  B 
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Do,  for  j=l  to  10 

(i)  Enter  energy  of  group 


j,  from  above  (b).  Press  R/S 
Note:  If  the  energ:  of  group  j  is  less  than 
12  keV,  then  repeat  this  step  for  the 
next  energy  group.  If  the  energy  is 
is  12  keV  or  higher,  then  continue 
with  (ii). 

(ii)  Enter  Kalansky's  coefficients  (Appen¬ 
dix  B)  for  the  nearest  energy  BELOW  the 
energy  of  the  group  (Aj , A2 ,Cj ,C2 ) , 
pressing  R/S  after  each  coefficient. 
Repeat  the  process  for  the  nearest 
energy  ABOVE  the  energy  group. 

(iii)  Enter  energy  corresponding  to  the  co¬ 
efficients  BELOW  the  group  energy, 
press  R/S.  Do  the  same  for  the  energy 
corresponding  to  the  HIGHER  energy. 

(iv)  Return  to  step(i)  unless  this  was 
group  10  in  which  case  go  to  three. 
Note:  Table  VI  identifies  the  data  registers  used  for  each 
coefficient  in  (ii)  above.  As  long  as  an  error  in  entering 
a  coefficient  is  detected  prior  to  entering  C2  high,  then  the 
correct  values  can  replace  the  incorrect  values  by  storing 
the  correct  value  in  the  applicable  storage  register. 
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TABLE  VI 

Data  Registers 

Coefficient 

Data  Register 

Ai  low 

24 

A2  low 

25 

Ci  low 

26 

C2  low 

27 

Ai  high 

28 

A2  high 

29 

Ci  high 

30 

C2  high 

31 

However,  once  C2  high  is  entered  for  each  group,  the  pro¬ 
gram  is  irrecoverable.  If  an  error  has  been  made,  the  pro¬ 
gram  must  be  re-executed  from  2,b,  above. 

3.  Enter  card  side  five.  Card  side  six  should  also  be 
read  in  if  either  the  4Trr2  fluence  or  fluence  cal¬ 
culations  are  desired.  Card  sides  five  and  six  are 
composed  of  the  program  on  pages  61  through  63, 

After  card  read-in 
Press  E 

Read  Transmission  Coefficient 
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TABLE  VII 


Data  Register  Entries 


Quantity 


Data  Register 


30402440 


13422200 

17310000 


Optional 


Press  E' 


Enter  source  yield  in  kilotons  Press  R/S 
Enter  x-ray  fraction  of  source  yield  Press  R/S 
Read  4irr2  fluence  in  calories 


Press  D" 

Enter  range  to  receiver  in  centimeters  Press  R/S 
Read  fluence  at  range  r  in  calories/cm 
These  instructions,  while  certainly  complicated,  are  hope¬ 
fully  clear  enough  to  enable  the  reader  to  use  the  program. 


Programming  and  Data  Registers 

Prior  to  entering  the  program  to  be  executed,  the  part 
tioning  of  the  TI-59  must  be  adjusted  to  639,39.  Addition¬ 
ally,  the  quantities  shown  in  Table  VII  should  be  stored  in 
their  respective  data  registers  prior  to  executing  the  pro¬ 
gram  or  recording  it  on  magnetic  cards.  The  other  26  data 
registers  are  used  as  working  registers. 


Appendix  G:  Program  PROMPT 


Program  PROMPT  is  used  to  calculate  the  prompt  radi¬ 
ation  fluences  or  doses.  The  applicable  coefficients  from 
Table  IV  must  be  input.  The  FORTRAN  PROMPT  listing  is 
shown  on  page  71.  A  TI-59  program  listing  is  shown  on 
pages  72  and  73 <  The  HTI  group  of  programs  does  not  con¬ 
tain  a  similar  program. 

FORTRAN  PROMPT 

The  input  variables  for  the  FORTRAN  version  are 
I  «  1  implies  a  thermonuclear  calculation 
2  implies  a  fission  calculation 
\  K  “  1  implies  a  neutron  calculation 

2  implies  a  secondary  gamma  ray  calculation 
A,B,C,D,E,F,G  are  taken  from  Table  IV 
XMI  is  the  mass  integral.  The  program  can  be  mod¬ 
ified  to  calculate  the  mass  integral. 

Note  that  the  output  statement  says  '4ttR2  FLUENCE'  ; 
however,  the  result  may  be  4ttR2  Dose  depending  on  the 
particular  problem  being  worked. 

TI-59  PROMPT 

To  use  the  TI-59  program  listed  on  pages  72  and  73,  15 
quantities  must  be  stored  prior  to  execution.  The  quanti¬ 
ties  and  their  storage  locations  are  shown  in  Table  VIII. 
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TABLE  VIII 

Prompt  Data  Registers 

Quantity 

Data  Register 

M.I. 

1 

A 

2 

B 

3 

C 

4 

D 

* 

5 

E 

6 

F 

7 

G 

8 

Slant  Range  (cm) 

9 

576357000 

11 

2127411731 

12 

1517003235 

13 

16323617 

> 

14 

212741 

15 

1731151700 

16 

3235001632 

17 

3617000000 

18 

Once  these  quantities  are  input,  pressing  B  will  cause 
the  program  to  execute.  As  in  the  FORTRAN  prompt  program, 
the  output  will  always  be  titled  4ttR2  FLUENCE;  however, 
the  answer  may  be  4nR2  DOSE  depending  on  the  input 
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coefficients. 


Output 

The  output  from  both  PROMPT  programs  is  in  the  form  of 
4nr2  fluence  or  dose  per  source  neutron.  Therefore,  to 
find  the  fluence  or  dose  at  some  point,  the  number  of  source 
neutrons  must  be  known. 


PROGRAM  PROMPT 
10  CONTINUE 

READ(#,#,END=999)  I » K  » A ,  B  >  C ,  D ,  E » F*»  G  >  XM I 
PRINT#,  '  ' 

PRINT#, 'THE  MASS  INTEGRAL  IS?  ',XMI 
I F ( I , EG . i ) THEN 

PRINT#, 'THERMONUCLEAR' 

ELSE  IF< I »EQ.2)  THEN 
PRINT#,  'FISSION' 

END  IF 

IF<K .EG. 1 >  THEN 

PRINT#, 'NEUTRON  CALCULATION' 

ELSE  IF (K . EQ . 2 )  THEN 

PRINT#, 'SECONDARY  GAMMA  CALCULATION' 

END  IF 

•  H  «  A  +  <B#XMI )  +  <C#XMI#XMI )  +  (D#<XMI##1 .5) >  + 
#(E#SGRT(XMI) )  +  <F#<XMI##( 1./3. > > )  +  ( G*ALOG ( XMI ) ) 
HH  «  EX  P  <  H  > 

PRINT#, 'THE  4  PI  R##2  FLUENCE  IS!  '  ,HH 
GO  TO  10 
999  STOP 
END 


fcHlS  PAGi,  u  QUALITY  FRACUGAjiLR‘ 
KROM  COrY  /UiuNISllEL  TO  DOG 
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Program  PROMPT 


000 
00  i 
002 
003 
004 
005 
006 
00? 
QOS 
009 
010 
oil 
012 
013 


LBL 
S' 


i  2 


0 

o 

o 

..I 

0 

4 

0 


uo 
00 
03 
00 
04 
00 
02 
04 
04 
00 

69  Op 
04  04 

014  43  ft CL 

Oi-5  01  01 

69  OP 
06  06 
69  UP 
00  00 


4 

4 

0 


43 


016 
01? 

01 S 
019 
020 

021  11  11 
022  69  DP 

01  *oi 
43  ft  CL 
12  12 
69  DP 
02  02 


023 

024 

025 

026 

02? 

023 

029 

030 

031 

032 

033 

034 


43 

13  13 

69  UP 
03  03 

43 

14 

69  DP 


14 


035 

036 

•03? 


033 

039 

040 

041 

042 

043 

044 

045 

046 

04? 

043 

049 

050 

051 

052 

053 

05W 


055 

056 

05? 

053 

059 

060 


04  04 

69  OP 
05  05 

6-9  OP 
00  oo 
53  ( 

53  e 
43 

02 

OP 
C»  •  J 

53 
43 
03 
65 
43 

01  01 
54  'i 
35  + 

53  < 

■*3  RCL 


02 

+ 

(’ 


r  ! 

03 


01  01 


33 

65 

43 


Vg 


RCL 


04  04 

54  > 


062 

063 


35 

53 


0?0 

071 

*  072 
0?3 
074 
0?5 
0?6 
07? 
078 
0?9 

•  080 
0  81 
n  o  o 

yoii 

083 
034 
035 
036 
03? 
033 
089 
09  0 
091 
092 
093 
094 


54 
65 
43 
05  0 

) 

+ 

y 

RCL 
01 
rx 

X 

RCL 
06 


54 
35 
53 
43 
.  01 
34 
65 
43 
06 
54 
85 
53 
53 


43 

01 
2  c; 


095 


4 
9.'*! 

03 

03 

03 

03 

03 

03 

54 


■b 
< 

< 

RCL 

01 

yx 


**v 

•v 

3 


064 

065 

43 

01 

RCL 

01 

098 

099 

65 

43 

•* 

RCL 

066 

c* 

^5 

y  v- 

100 

0? 

r 

06? 

01 

1 

:oi 

54 

V 

068 

93 

102 

85 

+ 

069 

05 

• 

5 

103 

CO 

V».« 

c 

V 

iri4  , 

•+  *  < 

RCL_ 

1 05 

1 06 
10? 

1 08 
109 

no 
in 
112 
X 13 

114 

1 15 

116 
11?  99 


01  01 
*  *'  O  (  *  *  *  * 


65 

43 

08 

54 

54 
•*' 
w  <L» 

’T* 
**•  V 

54 


RCL 

08 

5 


JNV 

LNX 


1  •  o 

i  l  O 

;  i  9 
120 
121 
122 
123 
. ;  24 
125 
?. 26 
12? 

1  :>  *:• 

4  £(>„i 

•"  **  *' 

1 30 

131 

132 

*  '*!»  '*i 

»•  o  »#i 

1 34 

1 35 


**3  STD 

1 0  1  o 
PRT 

RIiV 


OO 

?c* 


•:*  **' 

1  T  1 11  | 
*  It  V 

52 

EE 

43 

RCL 

15 

15 

69 

OP 

01 

01 

43 

RCL 

16 

*  «« 

69 

DP 

02 

02 

43 

RCL 

.!? 

1? 

69  I 

DP 

03 

03 

43  &}'■} 

18 

IS 

136 

4  *T« "? 

*  v* » 

138 

1 39 


69  or 
04  0 4 

69  DP 
05  05 

69  np 


&HIS  PAGE  IS  BEST  QUALITY  PiUCXLftASIi 
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I 


!  i 

I 


140 

00 

00 

141 

53 

< 

142 

53 

( 

143 

53 

144 

04 

4 

*+ 

145 

65 

x 

146 

83 

147 

65 

•.t 

148 

53 

< 

149 

43 

RCL 

150 

09 

09 

151 

33 

i  *  ■% 

152 

54 

> 

153 

54 

> 

154 

35 

1  /X 

155 

54 

> 

156 

65 

X 

157 

43 

RCL 

•158 

10 

10 

159 

54 

N 

160 

cm 

r  •* 

FR7 

161 

Q-? 

|U 

RTN 

162 

00 

0 

* 


S  K  mhm 

*  u;rx  ifUfu.ioiw.y  IODDO 
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